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■ INTRODUCTION 

o 
o 

(N . . . 

Statistical models in general apply to phenomena which appear in a sufficiently large number and can be repeated 
independently and (in principle) indefinitely. Whether high energy physics and in particular heavy-ion collisions 
provide such phenomena, is a question to be discussed on its own right. We shall try to contribute to an answer to 
[ — ■ this question here by reviewing some basic concepts of statistics, kinetic theory and thermodynamics, with the aim 
, to clarify the limits of the application of statistical models. 

, First we review the central limit theorem of statistics^, as the basic law governing the result of very many 

■ independent influences. Then kinetic theory, designed to describe the way towards equilibrium and its maintenance, 
, is discussed. By doing so we point out quite a few modern applications of statistical methods describing a stationary 

■ state, which obeys non-extensive thermodynamical rulesQ. A review of non-extensive rules, equilibrium one-particle 
, distributions and entropy density formulas follows. 

[ An underlying particular application we have in mind is a relativistic heavy ion collision!^. The first touch physics 
is probably dominated by nonlinear field dynamics and parton collisions. While some very energetic partons may 
escape in form of jets, in a central heavy-ion event most of the beam energy is transformed into a compression of 
' . • the nuclear matter and production of relativistic quarks and gluons. How a quark-gluon plasma in a thermal state is 
O , formed in these events and with which properties, is still an objective of the contemporary research. 
^ ' 

^ ' Assuming i) such a quark matter formation in an intermediate state, ii) ergodization of energy between many 
newly produced particles and iii) a relatively fast hadronization, the experimentally measured specific hadron spectra 
may reflect statistical, presumably even equilibrium thermal, properties of the precursor matter0|. Strong final state 
interaction on the other hand, a so called hadronic afterburner, may wash out spectral characteristics of earlier 
] stages of the evolution. Whether this is the case, can in principle be studied by checking quark coalescence rules or 
comparing hadron- and lepton-spectraj^ . Finally the late resonance decay during free streaming of hadrons changes 
the hadrochemical composition. Fortunately, some properties (e.g. the transverse momentum spectra) are influenced 
only partially (at their low end) by this. Since the relativistic energy is given by E = tot cosh ?/ with transverse 
mass Tot = -^/p^ + and rapidity y for a particle with mass to, the best way to study statistical equilibrium 
distribution of hadrons is the comparison of mT-spectra at rapidity y = for different particles. A universal behavior 
(TOT-scaling) indicates that the one-particle distributions depend on the energy only and not on all momentum 
components: a basic feature of generalized and conventional thermal distributions. 

In this lecture we review general basic ideas and mathematical formulas related to statistical models. By doing 
so emphasis is given to particular cases when the conventional picture of earlier textbook statistical physics does 
not apply, instead some generalized concepts have to be considered. Starting with the central limit theorem of 
statistics we continue with a discussion of the Brownian motion in the framework of the Langevin and Fokker-Planck 
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equations. Here the very level of generality is identified that leads to a Tsallis distribution instead of the Gibbs- 
Boltzmann oneQ. The potential for further generalization is pointed out, too. Then a discussion of the Boltzmann 
equation follows introducing to modern generalizations capable to describe non-exponential equilibrium one-particle 
distributions!^ llo|. We close with some remarks about how certain a thcrmodynamical state may be reconstructed 
from the observed one-particle energy spectra. 



STATISTICS: THE LAW OF BIG NUMBERS 



There is a mathematical property behind the applicability of statistical physics: many "normal" distributions of 
probability fold to a Gaussian with a width scaling down with increasing number of individual constituents. This 
property is expressed nicely in the central limit theorem^. An enlightening simple example of its action is given by 
the distribution of the sum of uniform random variables in a finite interval. Finally in this section an example not 
falling under the reign of the central limit theorem, the Lorcntzian distribution, is presented |lll . 

Let Xi be random distributed according to Wi{xi). We are interested in the distribution of a scaled sum of n such 
variables: 



/n / n \ 

i=l V k=l I 



(1) 



Here we assumed that the joint probability of having n values is a product of the individual probabilities; this is the 
requirement of statistical independency. With this assumption the Fourier transform of the seeked probability is a 
product of properly scaled Fourier transforms of individual probabilities: 

/n 
dxe^^'P^ix) = n*.(an^)- (2) 
i=i 

From the Taylor expansion of lnP(A:) around fc = one obtains the central moments (correlations). For the ^-th 
moment the scaling law, tri^' = a\jia^^^ applies with being the finite average of the given central moment of the 
individual distributions. Whenever af^ = for all i < £a values and af °) is finite, all the higher moments can be 
scaled down in the folded distribution by choosing a„ cx ri"^/^". The resulted scaling, an"^ = n^~^/^°a^^\ in the 
ri ^ oo limit leaves us with only one nonzero central moment, the £o-th one. 

Usually the £ ~ central moment is zero due to the normalization of the probability: InP(O) = In 1 = 0. The first 
moment = 1) can be made zero due to symmetry or by a simple shift in the variables Xi. The second central moment 
is then the first nonzero value. As a consequence the resulted distribution of the scaled sum has a second moment 
with all higher moments vanishing, therefore InP(fc) is quadratic in fc, so P{k) and with that P{x) is Gaussian. 

A nice, simple example is given for Xi-s uniformly distributed in the interval (—1,+!). The distribution of 
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(3) 
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approaches the Gaussian: 




lim P„(fc) = hm ( -^^"(^vg7^) j ^ cxp(-fcV2). (4) 



A counterexample is given by the Lorentzian distribution, having a Fourier transform Wi(fc) = exp(— |fc|). Now 
there is a non-Gaussian limiting distribution with an altered scaling for 
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Xi. (5) 
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It is itself a Lorentzian: 

lim P„(fc)= Km fe-l'^'/")" = cxp(-|fc|). (6) 

n — 'oo n — ^oo \ / 

This is a special case of the Levy distribution [l^. 



KINETIC THEORY 



In a brief review of concepts distilled from the kinetic theory approach to thermodynamics we shall rely on the notion 
of "noise" heavily. The sum of random influences (forces) is itself a random variable. It shows relationship to the sum 
of random numbers; its distribution under quite general circumstances can be considered as a Gaussian distribution. 
The independency of individual influences is assumed first of all in time: such a view deals with uncorrelated stochastic 
changes in the parameters describing a physical system. 

The simplest physical theory of such a system is that of the Brownian motion [l^. considering a free, massive particle 
under the influence of forces acting as an uncorrelated, Gaussian (white) noise. A description of the motion of such a 
particle is given by the classical Langevin equation, or equivalently a statistics over possible such motions is described 
by the Fokker- Planck equation. A balance between damping and accelerating forces leads to a stationary state, with 
vivid microscopical dynamics, but macroscopically (on the average over many particles) presenting a thermodynamical 
equilibrium state. General statements about this balance are comprised in the fluctuation-dissipation theorem. In 
the next section we review generalizations of the Boltzmann equation, a somewhat more complicated kinetic theory. 
The effects of such generalizations on the the concept and definition of entropy and on equilibrium distributions will 
be an issue of a further section. 



General Langevin problem 



Let us consider a simple, one degree of freedom motion. The change of momentum p in time is given by a force 
depending on this momentum and on a noise variable z: 

p = F{p,z). (7) 

Following the method pioneered by Ornstein and Uhlenbeck|l4l| a distribution of many possible values of p at a time t 
is considered. This /(p, t) distribution governs average values of a smooth but otherwise arbitrary test function R{p). 
The same integral over p can be expressed at the time t -\- dt assuming a statistical average over the noise z in the 
time interval passed since t: 



dpR{p)f{p,t + dt) ^ jdp{R{p + dtF{p,z)))f{p,t). (8) 

One assumes that the averaging of the force F over the noise z gives the result: 

(F) = -G(p), 

{FF)-{F){F) = 2D{p)/dt. (9) 

The above scaling of the correlation with dt follows from the Gaussian nature of the noise z. Expanding now the 
equationQ up to terms linear in dt one arrives at: 

dpR{p)^{p,t) = [dp[^G{p)R\p)+D{p)R"{p)]f{p,t). (10) 



dt 

After partial integration and considering arbitrary R{p) one gets the Fokker-Planck equation: 
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Particular Langevin problem 



The above Langevin and Fokker-Planck problem is still quite general. Damping and diffusion coefficients depend on 
the momentum p in a general way. Ergodization in phase space is achieved on the other hand when constant energy 
surfaces are covered. In such a situation the distribution /, as well as the coefficients G and D (the latter related to 
the noise), depends on the energy E{p) only. Note, however, that they are not constant. 

Such a particular Langevin equation is given bv[l^ 

p = z-G{E)^ (12) 

containing an energy dependent damping proportional to the general velocity v = dE/dp, and a zero- average noise, 
{z{t)) = 0; with a correlation 

{z{t)z{t')) = 2D{E)5{t - t'). (13) 

The Fokker-Planck equation contains in this case the factors D{p) — D{E) and G{p) = ~G{E)dE / dp. Its stationary 
solution is given by 

This result is not readily the Boltzmann-Gibbs distribution, cxp(— iJ/T), only in the case of energy-independent 
damping and noise coefficients. A general equilibrium is able to feature almost any other distribution of the energy 
of a single degree of freedom picked out of its environment. Instead of a constant temperature, T, in the general case 
a sliding inverse logarithmic slope is characteristic to such states. From \/T{E) — —d\nf{E)/dE its relation to the 
damping and diffusion coefficients follows: 

The low-energy limit of this expression, pretending as G{E) and D{E) were constant, leads to an experimentally 
feasible definition of the Gibbs temperature, Tcitbs = ^(0)/G(0). From the viewpoint of the Brownian motion 
another temperature may be used, the Einstein temperature TEmstein = ^^'raE^oo E>{E)/G{E). None of these two 
approximations arc, however, coincident with the sliding slope of particle spectra given by ca. H15|l . 

An important and historically mostly considered particular case is given by the constant slope distribution, the 
Boltzmann-Gibbs distribution: T{E) = T. A modern, non-classical distribution, the Tsallis distribution seems to be 
the next simplest, having a linear inverse slope ~ energy relation: 

T{E)^T/q + E{l-l/q). (16) 

The corresponding equilibrium distribution in this case turns to be an exponential of a logarithm, which is a power-law: 

1 / E\^ 
fip) = - l + {q-l)-] (17) 



Z \ 

It has the interesting property, that the parameter T is the fixed point of the sliding (linear) slope: T[T) = T. This 
parameter shall be referred to as the Tsallis temperature. 



Fluctuation-dissipation theorem 



In a realistic system there are many microscopic degrees of freedom to be considered. Denoting a point in the 
67V-dimensional phase space by pi, the Langevin equation can be implemented in the form 

K = (S^J^G^,)WJE + Zi. (18) 
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The symplectic coefficient, 5',^, does not cliange the total energy of the system, E{pi), it causes conservative motion 
inside a given energy shell only. For the sake of energy distribution it is therefore not interesting and shall be omitted 
in the followings. The damping and dissipation terms with the respective symmetric coefficients Gy and D^j keep 
balance on the long term. An ergodized equilibrium distribution can be a single function of the energy E only, 
therefore these coefficient matrices have to be connected by a single function of energy, too. This gives rise to a 
general fluctuation dissipation theorem: 

D,,iE) = T{E) {G.,,{E) + Dl^{E)) . (19) 

It is highly nontrivial that two high-dimensional matrix functions of the phase space coordinates pi would be related 
by a single scalar function of energy only! Recalling that T{E) is the inverse logarithmic slope of the equilibrium 
distribution, the fluctuation dissipation relation can be expressed using f{E), too: 

oo 

Du{E) ^ I G,,{x)f{x)dx. (20) 

E 

Still, quite general diffusion and damping coefficient matrices are allowed, but Dij is connected to Gij via an energy 
dependent scalar, the equilibrium energy distribution, f{E). Particular cases of this relation are i) the use of the Gibbs 
distribution, f{E) (x exp{—E/T), leading to Dij = TGij with constant matrices (the usual textbook case), or ii) the 
use of a Tsallis distribution (|17() giving rise to Dij{E) = (T + (g — l)E)Gij with constant Gij but linearly energy 
dependent diffusion coefficient, Dij. A further interesting case is presented by assuming that both the damping and 
the diffusion coefhcient matrix is a linear function of the energy, but their ratio (the Einstein temperature) is constant. 
This assumption is typical to field theory calculations. Due to ea. H19() the sliding slope is, however, not constant; it 
rather interpolates between a linear rise at low energy and a constant at high energy. From D = TG = ^T{1 + E/Ec) 
it follows 1/T{E) = 1/T +1/{E + Ec). 



Additive and multiplicative noise 



Another approach to reach a non-exponential stationary energy distribution from a kinetics described by the 
Langevin equation considers a stochastic damping coefficient. The Langcvin equation is kept linear, 

P^C-IP, (21) 

but, more general than in the classical approach, both ( and 7 are stochastic variables With constant mean 
values, (C) = F and (7) = G, and white-noise correlations, 

(7(t)7(t')>=2G5(t-t'), (C(t)7(t'))=2S<5(t-t'), {Cmt'))=2DS{t~t'), (22) 

the equivalent Fokker-Planck equation Hll|l contains a damping factor G{p) ~ Gp — F and a diffusion factor D{p) = 
Cp^ — 2Bp + D. For B the two noisy coefficients are cross-correlated. For a single degree of freedom the stationary 
distribution can be obtained analytically: 

'''' 

with the power v = 1 + G/2G, the exponent factor a = GB/C — F and the variable 9 = \J DG — S^. For F = a 
characteristic momentum scale in this distribution is given by p\ = D /G ^ the ratio of the additive and multiplicative 
noise widths. For F — Q and B — 0, i.e. for two independent noises, the Tsallis distribution arises as a stationary 
solution: 

fip) = A(^l + ^p'^ \ (24) 
Utilizing the energy formula for a free, massive, non-relativistic particle it reads as 

= + (9-1)1) (25) 
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with the Tsallis index q — 1 + 2C/G and the Tsallis temperature coinciding with that of the classical Brownian motion, 
T ~ D/mG. In the small momentum limit, p <^ pc the Tsallis distribution is nearly Gaussian, independent from 
the properties of the multiplicative noise, / = Acxp{—Gp^ /2D) (or expressed with the energy a Gibbs distribution, 
/ = Aexp{—E/T)). In the opposite hmit at high energy it is a power-law distribution, / = A{p/pc)~^^ = A{E/Ec)~'". 
It is interesting to note, that a definite relation arises between the energy scale, Ec, the temperature parameter, T, 
and the tail power, v: 

v=l + Ec/T. (26) 

This relation can be experimentally tested. 



NON-EXTENSIVE BOLTZMANN EQUATION 



The heart of kinetic theory is the classical Boltzmann equation. It docs not only describe dynamical evolution of 
large systems microscopically, resulting in stationary distributions on which thermodynamics can be based on, but 
it also offers a microscopical foundation to the key quantity entropy. It has, of course, also quite a few assumptions 
built in the theory; dropping one or other of them may lead to a generalization of the classical approach. The most 
recognized assumption, the micro-reversibility of the transition probability, establishes H-theorem and the definition 
of entropy; it should not be dropped. Less explicit assumptions, like taking the two-particle probability as a product 
of one-particle probabilities, or taking the total energy of a colliding pair as the sum of the respective one-particle 
energies for freely moving (asymptotic) particles, can be more readily generalized. 

We review these two generalizations of the Boltzmann equation in this section: the generalization of the product 
rule for probabilities (dropping statistical independency) leads to a non-linear Boltzmann equation (NLBeQ), while 
considering two- par ticle energies composed by an extended addition rule mounds in the non-extensive Boltzmann 
equation (NEBE[l(]j). Resulting stationary distributions, the H-theorem and the main characteristics of generalized 
thermodynamicscs following from this will be presented. 



NLBE: generalized product 



The general structure of the Boltzmann equation describes the evolution of the one-particle phase space density 
(interpreted as finding probability of a particle or a microstatc) via an integral over all possible transitions to and 
from other states: 

A = y W1234 (/3442 - /l2,34) ■ (27) 
234 

Here the dot denotes a total time derivative (Vlasov operator) comprising the essential evolution of the one-particle 
phase space density, /i. The indices 1234 refer to two particles before and after a microcoUision. The transition 
probability, W1234 reflects microreversibility and partner symmetry in the permutation of these indices. It is positive 
and it contains some conditions on conserving physical quantities; at least momentum and energy. 

^234 = M2234 6iipi +P2) - ip3 + Pi)) 5{Ei2 - E34), (28) 

with £'12 total two-particle energy before and E^^ after the collision. The particle density factors, /i2,34 and /34,i2 
weight the transition yields for a3 + 4^1 + 2 and for al + 2^3 + 4 process, respectively. 

In the NLBE approach the traditional simple additivity of energy is kept, 



E12 — El + E2, 



(29) 
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but the classical product formula of Boltzmann, /i2,34 = /1/2, or the supplement with blocking factors due to Uehling 
and Uhlenbeck /i2,34 = /i/2(l±/3)(l±/4) is generalized. The generalization still reflects particle separation property 
(particle 1 goes into particle 3 and particle 2 goes into particle 4 classically), but abandons the linearity: 

M34 = 7(/l,/3)-7(/2,/4). (30) 

This intends to simulate statistical correlations between initial and final states and a nonlinearity of transition yields. 
Further requirement is that the phase space density factor, 7 factorizes to a production factor a, to a blocking factor 
h (depending on the respective initial and final phase space densities only) and to a factor symmetric in both: 

7(x,y) = a{x) b{y) c{x,y) (31) 

with c{x,y) = c{y,x). 

The stationary state of the NLBE is governed by the ratio k{x) = a{x)/b{x). It is easy to see from the following 

derivation of the generalized H-theorem. We seek for a quantity called entropy in form of a one-particle phase space 

integral, S ~ J a{fi). The question is what cr(/) functional form guarantees macro-irrevcrsibility, i.e. S > 0. Using 

1 

the general Boltzmann equation H27|l one writes 

j iia'ifi) = y" wi234Ci3C24cri (03610462 - 01630264), (32) 

1 1234 

with the corresponding indices referring to arguments of the functions o, h and c. Now we explore particle permutation 
symmetries of this expression. After the exchange of particle 1 with 2 and simultaneously particle 3 with 4 we describe 
the same process. An exchange of the initial with the final state is done by 1 <-!■ 3 and 2 4 (micro-reversibility). It 
amounts to a relative minus sign in the total rate by exchanging gain and loss terms. Finally the combined operation 
of both also contributes with a minus sign in total. Dressing now the transition rate with factors carrying the same 
symmetry we use tl;i234 — ^123401302461626364/4. With the notation Ki ~ ai/bi for the production to blocking ratios 
we arrive at 

5= J W1234 (cr'i + 0-2 ^ CTs - cr4) (k3'«4 - K1K2) . (33) 
1234 

This already reminds to the structure of the H-theorem result of the classical Boltzmann equation. The correct 
entropy density can be read off as satisfying 

a'(/) = -ln«(/). (34) 

The stable equilibrium, where - due to ergodization - f{p) can be expressed as a solely function of the corresponding 
energy, satisfies a product rule for the K(/i)-s while the energy is additive. The only solution is that the canonical 
equilibrium of the NLBE is given by 

'^(/) = |exp(-i?/r). (35) 

As particular cases the well-known Boltzmann-Gibbs, Fermi-Dirac or Bose-Einstein distributions are recovered, but 
this result is far more general. The corresponding general entropy formula can be constructed founding a generalized 
thermodynamics with (in general) non-extensive entropy composition rules when merging large subsystems. 



NEBE: generalized sum 



Another, recently pursued way is to keep the statistical independency, /12.34 — /1/2, but to generalize the energy 
addition formula to a nontrivial composition rule 



Ei2 = h{Ei,E2). 



(36) 
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Rules not being a simple sum, h(x,y) ^ x + y, present a non-extensive energy composition. Latest in the thermo- 
dynamical limit, only associative rules are physical: h{h{x,y),z) = h{x, h{y, z)). Due to a mathematical theorem[r^ 
for functional equations the general solution of the associativity requirement is a strict monotonous mapping to the 
simple addition: 

X{h) = X{x) + Xiy). (37) 

This is unique up to a constant factor. 

As a consequence for any microcoUision X{Ei) +X{E2) = X{E3) +X{E4) holds. The stationary solution of NEBE 
is hence given by 

fip)^^expi-X{E)/T) (38) 

allowing again for a general, non-exponential energy dependence. It is noteworthy that to each associative non- 
extensive composition rule h{x,y) there exist a mapping function X{E). An in-medium dispersion relation, or 
quasi-energy is obtained this way: the total sum, 

X{Etot)^Y.^^^i^ (39) 

i 

is conserved by the NEBE. 

In order to connect the approaches NLBE and NEBE (and both with the non-extensive thermodynamics) the non- 
extensive composition rule of the energy has to be extended to composition rules of other, traditionally "extensive" 
quantities. Most important is the entropy composition rule. And, as putting a corner stone into the right place, the 
scaling law 

X{E)/T^Xs{E/T) (40) 

allows us to relate the equilibrium distributions by a one-variable functional relation, 

Zf,q{E) = eM-Xs{E/T)) = ^,-\cM-E/T)). (41) 

At the same time the generalized entropy density satisfies 

a'{}) = -\nK{Zf)^X;\-\nZf) (42) 

in equilibrium. This gives rise to an interpretation of the H-theorem for NEBE where the never-decreasing total 
entropy is given by the same strict monotonic back-mapping from the original Boltzmann-entropy: 

Sb = Xs{Stot)^ /"■^^''•^- ^^^^ 

The use of a non-extensive entropy composition rule, hs{x,y), mapped to additivity by Xs(t) leads to the following 
mapping of the individual entropy-density: s(/) = fX^^{— In /). The usual non-extensive entropy is then defined by 

ST = Jsif). 



THERMODYNAMICSES 



It is enlightening to list some particular cases of non-extensive composition laws and corresponding 
thermodvnamicscs|l8j. The trivial law, h{x,y) = x + y, is mapped by the identity, X{E) = E, and in (canonical) 
equilibrium the Boltzmann-Gibbs distribution, exp(— _B/r), emerges. The entropy density is given by s(/) = —/In/. 

Tsallis' non-extensive thermodynamics relies on the composition rule, h{x,y) = x + y + axy, connected to the 
mapping 



X,{S) = —\n{l + aTS), 



(44) 
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first presented by Abe[29j. The canonical equilibrium distribution, 

fe,{E) = ^il + aE)-'/''^ , (45) 

can easily be connected to the Tsallis distribution l|17|l by q = 1 — aT. Finally an entropy density following Tsallis' 
original suggestion (encountered also earlier by several authors0|): s(/) = "^f-p, can be obtained. 

An interesting endeavor is to consider the composition rule, h{x,y) = fx'' + y^)^^^ . Mapping the energy to its 
power, X{E) = (aE)^ /a, it gives rise in equilibrium to a Levy distribution|l2| . a.o. known from anomalous diffusion 
problems: 

f,g{E)^^exp{-iaE)'/aT). (46) 
The corresponding expression for the entropy is an incomplete gamma function. 

The pure multiplicative law, h{x, y) = axy^ which may be considered as the high-energy limit of Tsallis' construction 
with fixed energy scale E^, = 1/a, is mapped by X{E) = ^ \n{aE) and leads to a pure power-law behavior, Zf^q = 
{aE)~'^/°-'^ . The entropy formula due to Xs{t) ~ [l/aT) ln{aTt), s{f) = reminds to the original suggestion by 

Rcnvi[l9|. In fact defining 5*;^ = Xs{St) the (Tsallis) entropy is mapped to the extensive quantity: 

Sr = ^^ In / r- (47) 



Canonical and Extended Equilibrium 



Summarizing this part it seems to be useful to spell out the general conditions leading to the traditional and to the 
generalized treatment of canonical equilibrium. Three conditions are important: 

1. In the equilibrium state the energy shell is covered uniformly (ergodic assumption), 

2. The average energy exchange between the system under study and its environment is zero (not a driven system), 

3. The fluctuations in the system's energy are negligible compared to the average value of its energy (canonical 
limit). 

The third requirement has its physical basis in the fact that the characteristic range of forces acting in the system is 
much less then the system size. So in the thermodynamical limit of large systems, V ^ At., with a finite interaction 
range, ^, all surface effects (proportional to A) in a large volume, V , become negligible. In this case the entropy is quite 
generally additive, 6*12 = h{Si, S2) = Si + S2, and the equilibrium canonical distribution is of Boltzmann-Gibbs type, 
Zf = cxp{—E/T). Note that already the quantum statistical distributions violate this requirement: the exchange 
interaction (Pauli blocking or Bose enhancement) is long-ranged. 

In an extended equilibrium state, non-extensive thermodynamics and generalized kinetic descriptions deal with, 
the third requirement is not fulfilled. Either long range forces, or large fluctuations (related to each other) keep the 
non-extensivity parameter a ^ A£/V finite. Applying this idea to the time direction, also long-term memory effects 
spoil the traditional picture of canonical equilibrium and textbook thermodynamics. The unusual large fluctuations 
are treated in the framework of suDerstatistics[2(l| . where the traditionally intensive thermodynamical parameters, 
like temperature, have a probability distribution instead of a fixed value. This view, of course, also transforms the 
Boltzmann-Gibbs distribution into an arbitrary one: 

/e,(i?) = Jdprm^^e-^^. (48) 

For example, the Tsallis distribution can be obtained by considering a Gamma-distributed inverse temperature, 
(3 = 1/T. Some preview about a numerical simulation of NEBE for the Tsallis distribution can be obtained from a 
poster presented by G. Purcsel at this conference. 
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HADRON STATISTICS 



Statistical models have been often applied to hadron physics. Starting with Rolf Hagedorn's statistical model of 
meson resonances |2lj. several attempts occured to describe hadron multiplicities in elementary collisions by means of 
statistical distributions. The very idea of a phase transition between confined and deconfined quark matter relies on 
traditional equilibrium thermodynamics. With the dawn of relativistic heavy ion physics a search for the nuclear and 
quark matter equation of state began assuming a local thermal equilibrium in an otherwise exploding fireball. Such 
models and their predictions for experimentally observable properties will be reviewed in other lectures [23. Esf . 



Particle spectra and equation of state 



Here, following the above introduction to non-extensive thermodynamics, just some basic ideas about the inter- 
pretation of the power-law tails of nearly exponential spectra shall be discussed. Both exponential and power-law 
regions are observed in hadron spectra stemming from pp or heavy-ion collisions. The WT-scaling property, which 
reveals in particular after taking into account Dopplcr blue shift like corrections due to a transverse flow, indicates 
a thermodynamical origin. Although the traditional approach explains the power-law tail at very high pT values by 
pQCD calculations the non-extensive statistics provides a unified view for the whole spectrum. Certainly, jets also 
contribute to the hard part of light hadron spectra, but their angular distribution is far from uniform. In principle, 
these contributions may be filtered out. The question is, whether the rest still shows a (now statistical) power-law 
behavior. 

Another point can be made by inspecting the minimum bias pion pT-spectrum form RHIC AuAu collisions at 200 
GeV (Fig.l in Rcf. ^3)- To this a non-exponential fit can already be made at the p^-region between 1 and 4 GeV. 
The extrapolation of this fit almost coincides with the fit to the whole observed range between 1 and 12 GeV. So 
one concludes that the power-law behavior is not only a very hard scale physics. Furthermore in the non-extensive 
statistical approach there is a connection between the soft properties (temperature T) and the hard properties (power 
V, transition scale Ec) cf. eg . H2t)|) . 

The dilemma between the statistical and pQCD explanation can be well expressed in the following formula for the 
particle yield: 

oo 

^-^^ ^ J dcju;p{cj,k)fiuj/T). (49) 



The observed spectra, even not considering a collective flow or an extended source, are already convolutions of a 
spectral function, p{uj,k), and statistical-enviromental effects, /(cu/T). Both these factors may embrance non-trivial 
effects, distorting this way a naively expected Boltzmann-Gibbs exponential. Both may contain a further energy scale; 
the spectral function Kqcd, the thermodynamical weight Ec- Even in the case, when p would contain a very sharp peak 
at a given dispersion relation, Uk = w(fc), only, the result is a composite function, cxp(— w^/T), undistinguishable from 
a deformed thermodynamical equilibrium, exp(— X(|fc|)/T). Quite generally in this quasiparticle picture, whenever 
the effective dispersion relation, w(fc), contains parameters depending e.g. on the temperature, T, an assumption of 
an exponential equilibrium distribution like exp(— u;(fc, T)/T) questions the physical meaning of the temperature. If 
it is an enviromental parameter, the inverse slope is not necessarily equal to this, if it is an inverse slope, it is not 
necessarily the enviromental parameter governing an effective dispersion relation. 

Finally we note, that non-exponential distributions can also be interpreted in the framework of so called 
superstatistics|23. Here a distribution of the intensive thermodynamical parameters is assumed instead of a fixed 
value. The physical reason behind maybe the finiteness of the observed subsystem, as well as extraordinarily large 
fluctuations not scaling down with the system size alike the ones governed by the Gaussian distribution would do. 
The Tsallis distribution, for example, can be given by a continous distribution of the inverse temperature according 
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to a Gamma distribution: 



oo 



(l + 3;/c)-(^+i) = jdtt^e-'e--"^. (50) 



There are several pos sibilities to intrepret such a distribution. Among others a heat conduction equation with 
multiplicative noise |25l or taking into account an energy imbalance in two-body collisions due to a presence of further 
participant agents [23 leads to the required result. 



Limiting temperature with Tsallis distribution 



Another probe for the Tsallis distribution than the one-particle pT-spectra may be given by considering the spectrum 
of heavy hadronic resonances. An exponentially growing mass spectrum, originally proposed by Hagedorn and recently 
checked again latest experimental data in Ref.|23, with its famous consequence of having a limiting (or Hagedorn-) 
temperature for such a system, can be reconstructed on the basis of Tsallis distributed quark constituents. This 
approach|28j assumes that the Tsallis distribution of the quarks and antiquarks is folded into mesonic and baryonic 
distributions of the conserved total energy satisfying X{E) = ^iX{Ei). In general for an A^-fold convolution of 
massless constituents with d dimensional momenta one easily obtains that the average energy satisfies 

(^(^)) - ^E:b^- (51) 

This expression diverges first as the temperature reaches the limiting value, Th = i?c/rf starting from zero (for positive 
Ec = 1/a, i.e. for repulsive modifications of the extensive energy composition rule). This way Hagedorn hadrons 
emerge from Tsallis partons. This model also explains naturally, why the baryonic and mesonic mass spectrum seems 
to have a different rise: they contain different polynomial coefficients in front of three exponential factors. 



CONCLUSION 



In conclusion we have reviewed basic ideas of statistical physics and kinetic theory which may lead to a non-extensive 
thermodynamics. In particular the role of a multiplicative noise in the linear Fokker-Planck and Langevin approach 
has been emphasized. Generalizations of the Boltzmann equation towards non-factorizing yield factors or non-additive 
energy composition rules were then shown to lead again to non-extensive entropy definitions and thermodynamics. 
Finally these ideas have been related to hadron spectra observed in relativistic heavy ion collisions. We pointed out 
that not only the omnipresence and mx-scaling of power-law tails in particle spectra, but also a possible interpretation 
of the Hagedorn spectrum of heavy resonances may support the presence of a non-exponential equilibrium distribution 
in hot quark matter. 
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